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Abstract 

Recently, the 0(a) and 0(a s In a) radiative corrections to the orthopositronium lifetime have 
been presented in closed analytical form, in terms of basic irrational numbers that can be 
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evaluated numerically to arbitrary precision [Phys. Rev. Lett. 101, 193401 (2008)]. Here, we 

present the details of this calculation and reveal the nature of these new constants. We also list 

explicit transformation formulas for generalized polylogarithms of weight four, which may be 

OO ' useful for other applications 
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I. INTRODUCTION 



Positronium (Ps), the electron-positron bound state, was discovered experimentally 
in 1951 [1]. Since that time a lot of attention has been paid to the determination of 
its properties, including lifetime, decay modes, and spectroscopy The experimental and 
theoretical accuracies achieved by now being quite high, there is little doubt that quantum 
electrodynamics (QED) is the only interaction in this system. In fact, thanks to the 
smallness of the electron mass m relative to typical hadronic mass scales, its theoretical 
description is not plagued by strong-interaction uncertainties and its properties, such as 
decay widths and energy levels, can be calculated perturbatively in non-relativistic QED 
(NRQED) Q], as expansions in Sommerfeld's fine-structure constant a, with very high 
precision. 

Ps comes in two ground states, 1 5'o parapositronium (p-Ps) and 3 Si orthopositronium 
(o-Ps), which decay to two and three photons, respectively. Here we are concerned with 
the lifetime of o-Ps, which has been the subject of a vast number of experimental and 
theoretical investigations. Its first measurement [s| was performed later in the year 1951 
and agreed well with its lowest-order (LO) prediction of 1949 [4]. Its first precision 
measurement [sj, of 1968, had to wait 9 years for the first correct one-loop calculation [gj], 
which came two decades after the analogous calculation for p-Ps [7J being considerably 
simpler owing to the two-body final state. In the year 1987, the Ann Arbor group |8J 
published a measurement that exceeded the best theoretical prediction available then 
by more than 8 experimental standard deviations. This so-called o-Ps lifetime puzzle 
triggered an avalanche of both experimental and theoretical activities, which eventually 
resulted in what now appears to be the resolution of this puzzle. In fact, the 2003 
measurements at Ann Arbor [9( and Tokyo [10], 

T(Ann Arbor) = 7.0404(10 stat.)(8 syst.) (is' 1 , 
r(Tokyo) = 7.0396(12 stat.)(ll syst.) (is 
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agree mutually and with the present theoretical prediction, 

r(theory) = 7.039979(11) /is" 1 . (2) 

The latter is evaluated from 

2 „/a\ 2 3a 3 



T (theory) = T 

where 



1 + A-H \na + B \n 2 a + C — In a 

7T 3 V7T/ 2fl 7T 



(3) 



r = -(7r 2 -9) — (4) 

y % 
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is the LO result. The leadin g lo garithmical 
terms were found in Refs. [Ill, 12 1 and Ref. 



y enhanced 0(a 2 lna) and (D(a 3 In 2 a) 



13| . respectively. The coefficients A 



-10.286606(10) 0, [ll|, Q, EI B = 45.06(26) H, and C = -5.51702455(23) [l7| 
were evaluated numerically in a series of papers. Comprehensive reviews of the experi- 
mental and theoretical status of Ps may be found in Refs. 18|, [l9[ . 

We note in passing that high-precision tests make Ps also a useful probe of new physics 
beyond he standard model. At present, there is strong interest in models with extra 
dimensions 20], which may provide a solution of the gauge hierarchy problem j2l| (see 
Ref. [22! for a review) . Some time ago, a peculiar feature of matter in brane world 
was observed in Ref. 23j, where it was shown that massive particles initially located on 
our brane may leave the brane and disappear into extra dimensions. The experimental 
signature of this effect is the disappearance of a particle from our world, i.e. its invisible 
decay. The case of the electromagnetic field propagating in the Randall-Sundrum type of 
metric in the presence of extra compact dimensions 24, E^J was considered in Ref. 25]. 
where it was shown that the transition rate of a virtual photon into extra dimensions is 
non-zero. This effect could result in the disappearance of a neutral system. In the case 
of o-Ps, such estimations for the invisible decay branching fraction 5(o-Ps — > invisible) 



19, 



261 ] range just one order of magnitude below the presently best experimental upper 
bound of 4.3 x 10~ 7 at 90% confidence level established by Badertscher et al. 27j]. Thus, 
this decay is of great interest for the possible observation of effects due to extra dimensions. 

In order to reduce the theoretical uncertainty in the o-Ps total decay width T(theory), 
it is indispensable to increase the precision in the coefficients A, B, and C in Eq. ([3]). This 
is most efficiently done by avoiding numerical integrations altogether, i.e. by establishing 
the analytic forms of these coefficients. The case of B is beyond the scope of presently 
available technology, since it involves two-loop five-point functions to be integrated over 
the three-particle phase space. In the following, we thus concentrate on A and C. The 
quest for an analytic expression for A has a long history. About 25 years ago, some of 
the simpler contributions to A, due to self-energy and outer and inner vertex corrections, 
were obtained analytically 28J], but further progress then soon came to a grinding halt. 
In our recent Letter [2^], this task was completed for A as a whole. The purpose of the 
present paper is to explain the most important technical details of this calculation and to 
collect mathematical identities that may be useful for similar calculations. 

An analytic expression for C is then simply obtained from that for A through the 
relationship [l3] 

A 229 

3 30 V ; 

which may be understood qualitatively by observing that the 0(a 3 In a) correction in 
Eq. ([3]) receives a contribution from the interference of the relativistic 0(a) term from 
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the hard scale with non-relativistic 0(a 2 In a) terms from softer scales. 

The structure of this paper is as follows. Section [III contains the well-known integral 
representation of the o-Ps total decay width as given in Ref. In Sec. IHH we show 

how to transform the contributing integrals to forms appropriate for analytic evaluation, 
which is carried out for the most complicated integrals, which are plagued by singularities, 
in Sec. IIVI More examples are studied in Sec. |Vj The final results for the coefficients A 
and C are presented in Sec. IVII Section IVIII contains a summary. In Appendix [A], we 
present the analytic results for all parts of the integral representation given in Sec. [TT1 
Appendix [B] contains useful representations of the ip function and the expansion of the T 
function about half-integer- valued arguments. In Appendix [Cj transformation formulas 
for generalized polylogarithms of weight four with different arguments are collected. 



II. DEFINITIONS AND NOTATIONS 




FIG. 1: Feynman diagrams contributing to the total decay width of o-Ps at 0(a). Self-energy 
diagrams are not shown. Dashed and solid lines represent photons and electrons, respectively. 

The 0(a) contribution in Eq. (JS]), Ti = T Aa/ii, is due to the Feynman diagrams 
where a virtual photon is attached in all possible ways to the tree-level diagrams, with 
three real photons linked to an open electron line, and the electron box diagrams with an 
e + er annihilation vertex connected to one of the photons being virtual (see Fig. [I]). Taking 
the interference with the tree-level diagrams, imposing e + e~ threshold kinematics, and 
performing the loop and angular integrations, one obtains the two-dimensional integral 
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representation [16] 



1 

ma 7 f dxi dx 2 dx 3 

r i = / 6(2 -Xi -x 2 -x 3 ) F(xi,x 3 ) +perm. 

36rr J xi X2 x 3 



(6) 



where x^, with < Xi < 1, is the energy of photon i in the o-Ps rest frame normalized to 
its maximum value, the delta function ensures energy conservation, and "perm." stands 
for the other five permutations of xi,x 2 ,X3. The function F(xi,x 3 ) is given by 

5 7 
F(xi,x 3 ) = g {xi,x 3 ) + y]gi(xi,x 3 )h i (x 1 ) + ^2g i (x 1 ,x 3 )h i (x 1 ,x 3 ), (7) 

i=l i=6 



where gi(xi, x 3 ) are ratios of polynomials, which are listed in Eqs. (A5a)-(A5h) of Ref. 
and 



16 



h^xi) = \n(2xi), h 2 (xt) = x f^9 1 , h 3 ( Xl ) = -^[C(2) - Li 2 (l - 2x0], 

»' xi 2x\ 



h*(xi) = ^"[3C(2) - 29% h 5 ( Xl ) = ^-61 (8) 

h 6 (xi,x 3 ) = = _ [Li 2 (rj[,0i) - Li 2 (r^,0i)l , (9) 
VX1X1X3X3 

M^i^3) = T7 = _ [2Li 2 (r+^)-2Li 2 (r-,^)-Li 2 (r+,0) + Li 2 (r-,0)], (10) 
where x,- — 1 — x,- and 



fl 4. /^l S 4. /^l /^l^S /^1^3 

#i = arctan W — , 9\ — arctan A / — , Pa = \ — — , Pb — \ > 

V x\ y xi v ^1^3 V ^1^3 



^(l±p A ), r± = ^7(l±^), r% = -2=. (11) 

sJX\ 



Here, C(2) = tt 2 /6 and 



1 

1 /* di 

Li 2 (r, 0) = -- / — ln(l - 2rt cos 9 + r 2 t 2 ) (12) 



is the real part of the dilogarithm [see line below Eq. (132 jl ] of complex argument z = re 



id 



301 ] . Since we are dealing here with a single-scale problem, Eq. ([6]) yields just a number. 

Although Bose symmetry is manifest in Eq. ([6]), its evaluation is complicated by the fact 
that, for a given order of integration, individual permutations yield divergent integrals, 
which have to cancel in their combination. In order to avoid such a proliferation of terms, 
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we introduce an infinitesimal regularization parameter 5 in such a way that the symmetry 
under for any pair i ^ j is retained. In this way, Eq. (JS]) collapses to 

1-5 l-<5 

T 1 = —— / dxi / —F(xi,x 3 ), (13) 

67H y J x x x 2 x 3 

where £3 = 2 — xi — x 2 - Note that we may now exploit the freedom to choose any pair of 
variables Xi and Xj (i ^ j) as the arguments of F and as the integration variables. 



III. INTEGRAL REPRESENTATIONS OF DILOGARITHMIC FUNCTIONS 

Obviously, the functions h^{xi : x 3 ) and h?(xi,X3) in Eqs. (jHJ) and ( |T0|) . respectively, 
give the most complicated contributions to IV In order to perform integrations involving 
these terms, it is useful to apply the integral representation of Eq. ffl~2]) to Li 2 (r^, 6 1 !), 
Li 2 (r^,^i), and Li 2 (r^,0). Let us first consider Li 2 (r^, #1). We see from Eq. (TTTT) that 
cos 6*i = y/x[ and thus 

l+PB 

Li 2 (r+e!) = -i j ^ln[l-x 1 t 1 (2-t 1 )], (14) 


where t\ = (1 + pfl)t. Then, the term D\ = Li 2 (rJ, 6\) — Li 2 (r^,6 l i) on the r.h.s. of 
Eq. (|T0l) . after the change t 2 = ti — 1, can be rewritten as 



Ps 

1 r dt 



D 1 = -- J ^-lnfl-^a -*!)]. (15) 



Finally, substituting t 2 = psVi, we obtain 



1 



1 /* c^i 

Di = --Vx7x\xsx^ —=- [In xi - In x 3 + ln(x 3 + x 3 t)] . (16) 

^ j Vt{xix 3 - xix 3 t) 

The residual term on the r.h.s. of Eq. (|T0l) . .D 2 = Li 2 (rJ, 0)— Li 2 (r^, 0), can be transformed 
in the same way yielding 

1 

D 2 = --V^l^i^3 / -7=7 [ln(gig 3 ) - ln(xix 3 ) + hit]. (17) 
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We thus obtain the following integral representation for h 7 (xi,x 3 ) [48]: 



/i7(xi,x 3 ) 



dt 



4 J y/t (xiXs — X\X 3 t) 



In ^1 _|_ 2 ln(x 3 + x 3 t) - In t 
x 3 x 3 



Exploiting the x\ <-» £3 symmetry of the coefficient g 7 (xi, £3) multiplying h 7 (xi, x 3 ) 
Eq. (|T8|) can be effectively replaced by 



h 7 (xx,x 3 ) 



dt 



4 y -y/t (0:1X3 — XiX3t) 



[21n(x 3 + x 3 t) - hit]. 



(19) 



Next, this expression, multiplied by g 7 (x\,x 3 ), is to be integrated over x\, x 3 , and t. 
Observing that the logarithmic terms in Eq. (fT9j) are independent of x%, we first integrate 
over x\ (for a similar approach, see Ref. [3l|). In order to avoid the appearance of 
complicated functions in the intermediate results, the integration over t in Eq. ( |T9i) is 
performed last. 

Using the same technique, we obtain the following representation for the function 

h & {x x ,x 3 ): 

1 



h 6 (xi,x 3 ) 



dt 



\/t(xiX 3 — X\X 3 t) 



[lnxi - In x 3 + ln(x 3 + x 3 t)], 



(20) 



in which the part proportional to lnxi and the complementary one are first integrated 
over x 3 and x\, respectively. The t integration is again performed last. 

In Sees. [IV] and El we discuss in more details how these integrations can be performed. 



IV. EVALUATION OF CONTRIBUTIONS WITH h 6 AND h 7 

We now discuss the evaluation of the most complicated integrals, namely those in- 
volving the functions h G (xi,x 3 ) and h 7 (x\, x 3 ). We denote the corresponding integrated 
expressions as I e and I 7 , respectively. They are both singular for 5 — > 0, so that the 
regularization of Eq. (TTHT) is indispensable. 

Let us first consider the contribution of the coefficient g 7 (xi,x 3 ) without the function 
h 7 (x\,x 3 ). It can be decomposed into two parts, as 

~ / \ 97\JEli%3) -sine/ \ . ~reg/ \ /01 \ 

g 7 (x l7 x 3 ) = 77; r = g 7 s (xi, x 3 ) + g 7 8 (x x , x 3 ), (21) 

XiX 3 {2 — Xi — x 3 ) 

where 

gt ns (x 1: x 3 ) = l X3{1 ~ X3) (22) 
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gives rise to the singularity upon integration over x\ and X3, while 

4 



g™ s ( Xl ,x 3 ) = — -3 + 9x 3 + ( 10)x!+ ( 

X3 \ x 3 / V 



2 - x 3 x 3 



— + 10 + 2x 3 — 



1 



+ 



52 



12 

— + 66 - 44x 3 + Uxl 



2 — xi — X3 



(23) 



2 - x 3 x 3 

remains finite, so that the limit 5^0 can be taken. A similar decomposition can be made 
also for ge(x\, X3). Specifically, performing the integrations over xi and X3 and taking the 
limit 5 —>■ 0, we have 

1-5 1-5 

6 / dxi 



dx 3 g? ns (xi,x 3 ) = 3 In 5+ - + 0(8), 



25 l-xj_+5 
1-5 1-5 



6 / dx\ 



dx 3 g r 7 eg (x!,x 3 ) 



1240 



264C(2) + 0(5). 



(24) 



26 



l-xj+5 

Observing that the presence of the functions h${x\, x 3 ) and hi{xx, x 3 ) does not change 
the singularity structure of the integrals over the variables £3, xi, and t in this order, the 
decomposition of Eq. (l2~il leads to 

1-6 1-6 



rsing ,-reg 
1 i "r 1 i ' 



r smg, reg 



6 / dx 3 



dx\ g 



~smg, rcg 



Xl,X 3 )/li(xi,X 3 ) 



(25) 



26 



l-a; 3 +5 



with i = 6,7. 

Our evaluation yields 

7-sing 



91n5 + 45 + -C 2 -yCs + 0(5), 



rreg 
J 6 



422 



+ C 2 ( 



1877 



1590/ 2 - 288Z? ) + 
35 



2 16 



r 



sing 



rreg 
1 7 



- 576Li 4 (-1 +—G 3 + 0(5), 
-91n5-36-yC 2 + yC3 + ^), 

297 + C2 ( - 222 + 486/ 2 J - —(3 + —(4 + — G 3 + 0(6), 



where 49] 



16 



21 



y/2 



(26) 

(27) 

(28) 
(29) 



G 3 = I2C2/2 - It - 3% 2 h. ~ W 2 l R + l 6 R - — Cs + 48 Li 3 -= 



+ 3 Re 



Li 3 



1 - y/2 



Li 3 



4 

1 + v 7 ^ 



V2 



(30) 
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As can be seen from Eqs. ( )26l) and (EE}, In 5 cancels in the sum I§ + I7. Here and in the 
following, we use the short-hand notations 



/ 2 = ln2, / 3 = hi3, l R = \n(l + y/2) 



Furthermore, 



Sn,p(>^) 



(n — ly.pl 



ln p (l - tx) In"" 1 t 



(31) 



[32) 



is the generalized Nielsen polylogarithm, Li n (x) = S n -i,i(^) the polvlogarithm of order 
n, and ( n = ((n) = Li n (l), with ((x) being Riemann's zeta function 30l. 132]. 

The result of Eq. (130]) , which is the most complicated part arising from the terms with 
i = 6 and 7 in Eq. ([7]), assumes a rather simple form when written as an infinite series, 



3 3 



14C 3 -24C 2 / 2 -iy 3^4" 

S3 S2 2 2 ^ r 2 , 

n=l v ' 



n + 2 



71+1 



(33) 



where ip( m \n) is the (m + l)-th logarithmic derivative of the T function, T(x) 
J °° dte - *^ -1 . We can now apply the well-known relations for T and if) functions, 



iP> 



n + 2 



-if)' 



T 2 (n) 1 
2r(2n) ~ ft) 
n + 1 



n 



where 



S± m {n) = ^2 



1 »4 



r±iv 



--C2 - S_ 2 (n) 



r 



(34) 
(35) 

(36) 



is the harmonic sum. Using Eqs. (13~4"1) and f[3"5"j) . the constant G 3 is rewritten in terms of 
so-called inverse central binomial sums, i. e. sums of the form 



~ii 

n=l \ n ) 



</>(ri), 



where 4>(n) is some combination of harmonic sums and factors 
number. Sums of such type were studied in great detail in Refs. 



(37) 



ike 1/n, and z is some 



33, 



34. 



35 



36 



It is known that, for the series in Eq. (j3"7|) . there exists a nonlinear transformation, 

(38) 



^/z - 4 - v/I 



v / ? 3 4 + yfz 
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which leads to great simplifications in many cases. The series in the new variable y 
does not have a binomial coefficient and can be summed, yielding expressions involving 
generalized polylogarithms S n , P (y)- 

Now we can explain the appearance of the prefactor 1 / y2 in front of G% in Eqs. (1271) 
and (129]) . Such a prefactor has not appeared in single-scale calculations so far. The point 
is that all inverse binomial series involving products of the factor 1/n and some function 
f(n) that is a combination of the ip function and its derivatives have the form (see, for 



example, Ref. [361 ]) 



where F(y) is some combination of generalized polylogarithms and y is defined by Eq. 
Note that Eq. (135]) contains the binomial sum S_ 2 (ra), which is related to the basic one, 
S-2( n ~ I); via 

S- 2 (n) = S- 2 (n-l) + ^f-. (40) 

Thus, the last term on the r.h.s. of Eq. (1401) leads to z = 4 in Eq. (1331) . which translates 
to y = — 1 via Eq. (!38l) . This term then cancels the term 14£ 3 — 24£ 2 / 2 on the r.h.s. of 
Eq. (133]) . For the term — £ 2 /2 — S-2(n — 1) on the r.h.s. of Eq. (135|) . we have z = —4 so 
that the variable y from Eq. ( 1381) assumes the value 

r = ^Zl. (41) 

This explains the appearance of the factor l/v^2 in Eqs. (127|) and (1291) . since (1 — y)/(l + 

y) = 1/V2. 

Finally, we can rewrite Eq. ( |30l) as 

G 3 = 21C 2 / r - ^ r 3 - 5/ r Li 2 (r) + 5 Li 3 (r) - 50 Si, 2 (r) + 4 Si, 2 (r 2 ) + 34C 3 , (42) 
where L = lnr. 



V. EVALUATING INTEGRALS FROM SERIES 



Let us now consider several typical integrals that arise upon the first integration 50]. 
Our first example of the remaining two-fold integrals reads 

i i 

r Ht r fir 

I ± = — ^ln[l T 4t(l-t)(l-x)]ln(l-x). (43) 
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Direct integration over t or x would lead to rather complicated functions in the remaining 
variable. Instead, we Taylor expand the first logarithm using ln(l — q) = — Y^n=i 1 n l n to 
obtain 



f^/^d-Op/fd-xj-bKi-*). 



(44) 



Now, the two integrals are separated and can be solved in terms of Euler's T function. 
Using 



J ^(1 - x) n ln(l -x) = -i>\n + 1), 



(45) 



we finally have 

n=l v ' n=l \ n / 

Clearly, in the cases of I + and J_, the argument z in Eq. ( 1371) is equal to 4 and —4, 
respectively. 

The case of I + is simpler and leads to a smaller number of constants. Indeed, we can 
use the results of Refs. 35|, |36j to obtain 



1 4 r 



n=l V n / 



E 

n=l 



2n^ ^2 



1 4 r 



3(2 



00 1 4^ 



15 



( 2 D 



n' 



n=l \ n J 

± + 8 Li 4 - 
3 \2 



4- 



(47) 



and so on. According to Ref. 36], after transformation to the variable y of Eq. (I38|) . 
we arrive at polylogarithms of argument —1, which are expressed in terms of alternating 
and non-alternating Euler-Zagier sums, such as C(i fl ) = S^li( = ' = l) r V na > C(=t a ) = '=^) — 

E™=i E^° =m+ i(±l) n (±l) m /K^), etc. 

Let us now turn to the case of I_ in Eq. (1461) . The argument z = —4 gives y = r and 
leads to a new type of constants. Again using formulas from Ref. 36J, we have 



n=l V n I 



2n ^ n 2 



2 r ' 



00 1 f— zH n 1 

n=l V n / 



E 



1 (-4)' 



i ( 2 D 



-^ r 3 - + 4S 2 , 2 (r) - 4Li 4 (r) + 4/ 2 [Li 3 (r) - Cs] 

+ 4/ r [Li 3 (r) - Si, 2 (r) -/ 2 Li 2 (r)] -/ 2 [2Li 2 (r) + / 2 ]. (48) 
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With the help of the relations listed in Appendix EJ I± can be alternatively expressed as 

li 17. /T 



I + = -4C 2 / 2 2 - | + yC4 - 8Li 4 (J 1 , 
I- =(4- 7^2 + 2^2(2 + 5ljl 2 R - —l 2 R ( 2 - -^l R 
- M R Re 



4 Re 



(49) 




(50) 



It has been observed empirically that, at weight four, the terms that are not expressed 
through the usual Riemann zeta function ((a) often come in the combination 6 4 = l 2 (ll/3— 
2(2) + 8Li 4 (l/2) introduced by Broadhurst in Ref. |38l| . Examples include the three-loop 
QCD correction to the electroweak p parameter [39], the electron anomalous magnetic 
moment at three loops 40] , critical exponents in high orders of perturbation theory 4l) , 
;he heavy-quark contribution to the vacuum polarization function at four loops in QCD 
421 ] , and the matching conditions for the strong-coupling constant at five loops 43J . Our 
result for I + in Eq. (T4"9l exhibits a violation of this empirical observation. In fact, the 
non-zeta terms form some different combination there. 

Another class of typical integrals yields sums involving ip functions of half-integer 
arguments (see Appendix B), e.g. 



J4 



1 dt j d ln[l=F4t(l-t)(l-x)]ln(l-x) 

T J X x -2 





E 



(±4) n T 2 (n 



it 

E 



' 8n T(2n 
1 (T4) n 



, fn + 2 



2n\ n 2 



t in) 



-C 2 -S- 2 {n-l) 



n + 1 



-1^ 



rr 



(51) 



Following a similar strategy as above and using formulas from Sec. IIVt we may express 
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J± in terms of known irrational constants, as 



-9Re Li4(^j+Li 4 fAt^, +4 



Li4 i^ +Li4 fi±.^' 



+ M 3/2 + 5(2^2 + ^3 -16 Li 3 



+ 5 Re 



Li 1 I - — I + Li 4 



L13 



.2 — . 

Li 4 + Li 4 




(52) 



where G4, expressed with the help of the variable r defined in Eq. (1411) . is given in Eq. flA3j) . 

These results again contain various contributions of polylogarithms with argument 
y = —1, arising from terms of the form (— l) n /n 2 on the r.h.s. of Eq. (1511) for J + and 
terms of the form — C2/2 — S^n — 1) on the r.h.s. of Eq. (ISTT) for J_, and with argument 
y = r, arising from the residual terms. 

Unfortunately, not all integrals can be computed so straightforwardly. In more compli- 
cated cases, the integrations are not separated after expansion to infinite series. We then 
rely on the PSLQ algorithm 44j, which allows one to reconstruct the representation of a 
numerical result known to very high precision in terms of a linear combination of a set 
of irrational constants with rational coefficients, if that set is known beforehand. The ex- 
perience gained with the explicit solution of the simpler integrals helps us to exhaust the 
relevant sets. In order for the PSLQ algorithm to work in our applications, the numerical 
values of the integrals must be known up to typically 150 decimal figures. However, for 
some integrals more accurate determinations are required. The success of the application 
of the PSLQ algorithm also relies on the fact that only certain combinations of polyloga- 
rithms, like G3 in Eqs. f l30|) and f|42|) . G 4 in Eq. (IA3I) . and G4 in Eq. flA2j) are incorporated 
as independent structures. 



VI. RESULTS 



Finally, to get rid of complex polylogarithms, such as Li 4 [(l + y/2)/2], in the above 
formulas, we transform all polylogarithms to arguments of value below unity. To this end, 
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we need transformation formulas through weight four. Some of these formulas are listed 
in Appendix O After a laborious calculation, we obtain the final result for the one-loop 
correction 



2, 2 x , 56 19, 



901 2701, 253, 2 

k H h 

216 108 24 2 



11449 



432 



-(3 



_59_ 4 

576 2 



12983 
192 



251 



C 4 + — Li 4 ( - ) + G 4 + -G 



6V2 



(53) 



where the constants G 3 , G 4 , and G 4 are specified in Eqs. (1421) . (1A3I) . and (|A2|) . respec- 
tively. Transforming the polylogarithmic functions by means of the formulas given in 
Appendix [Cj we arrive at the form of Ref . [3] , 



2, 2 v , 56 19, 
-(tt 2 -9)A= — + — / 



901 2701 , , 11449 _ 253, ,2 913 , ,2 251, 4 

C2 C2/2 H C3 H (2^2 H C2J 2 H ^2 

216 s 108 S 432 S3 24 S2 2 64 S2 3 144 2 



27 6 

83 4 91 11303 

256 3 6 192 s 

21 49 2 7 3 35 4 35 581 fl 

4 s 16 s r 16 r 384 r 8 s 16 S \3 

- ^-l 2 Li 3 (-r) - ~l r Li 3 (-r) + S 2 Li 3 (r) + S r Li 3 (r) 

Z Z 4 o 



249 
132" 
63 



Li 4 - 



249 T . 

H Li 4 

16 



251 _ . 

+ ^r Ll4 

6 



+ 7Li 4 (-r) -7S 2 , 2 (-r) 



T . / % 63 . . 7 
r L l4 (r) + T S2,2(r) + - 7 = 



7 15 

-C2^r lr h Li 2 (r) 

2 72 r 6 



+ - Li 3 r - ~r Si,2 (r) + - Si, 2 (r 2 ) + — Cs 
6 3 3 3 



(54) 



where r is given in Eq. (Hlj) . 

From Eqs. (1541) and (15"]) , A and C can be numerically evaluated with arbitrary precision, 



A = -10.286 614 808 628 262 240 150 169 210 991 253 179 644 007 490 228 232 410 . 
C = 5.517 027 491 729 858 271 378 866 098 665 005 181 944 001 421 860 702 103 921 



•(55) 



These numbers agree with the best existing numerical evaluations ly, [l7j] within the 
quoted errors. 



VII. CONCLUSION 

We presented the details of our evaluation of the 0(a) and 0(a 3 In a) corrections to the 
total decay width of o-Ps, i.e. of the coefficients A and C in Eq. d3J), respectively, which 
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had been presented in our previous Letter [29] in closed analytic form. We discussed the 
nature and the origin of new irrational constants that appear in the final results. They 
were shown to be some particular cases of inverse central binomial sums and corresponding 
generalized polylogarithms. These constants enlarge the class of the known constants in 
single-scale problems. 

The 0(a 2 ) correction B in Eq. ([3]) still remains analytically unknown. 
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APPENDIX A: DETAILED RESULTS 



In this appendix, we present separate results for the integrals of Eq. (1251) with 
i = 0, .... 7. Note, that not all of them a finite in the limit 5 — > 0. We have: 



J = 204 - 142C 2 , 

362 1273 
I, = 51 + 90Z 2 - 228C 2 + — ( 2 k + -^Cs, 



h 
h 



-40 - - 72C 2 / 2 + 42C 3 - ^=G 3 , 

5 V2 



144( 2 In 5 - 59 + 24/ 2 + ( 2 



— + 371/ 2 - 294lt) + 52C 3 - 17Z| - 



16 



408 Li 4 ( - ] + 36G 4 , 



328 



+ C2 ( - 252/ 2 + — 1' 2 ) + 35C 3 4 



380 



+ 1026Li 4 ( - ] + 27G 4 , 



783. 



279 
16 



1863 



109 



•144C 2 In 5 - 120 + 6C2 (-3 + 6/ 2 + 95Z|) - 357C 3 + — 1\ - 1398C 4 



+ 1464 Li 4 



287 

91n5- — + C2 



+ 36G 4 , 
/3781 



6 



576 Li 4 



35 

7! 



2559 

- 1590/ 2 - 288/^ ) + 1328Cs - 24/^ + (4 



G, 



h = -91n<5 + 261 + C 2 



471 



315 



24 



+ 486/ 2 - 252C 3 + — C 4 + -=G. 



16 



V2 



where G3 is given in Eqs. ( 1301) and (H21), and 



G 1 



913 2 581 . 

^3 + — C 2 Ll2 



n 249 
+ 12 



119 
~12 



Ga 



i^ + i^ 3 + ^ 4+5C4+C2 



9 



2Li 4 



Li 4 



7 fr ) + C3 



-3/ 



83 
256 



5 Z 

2-^ 



+ [9k + -Z r J Li 3 (r) + (-6/ 2 - 2/ r ) Li 3 (-r) - 9 [Li 4 (r) - S 2j2 (r)] 
+ 4[Li 4 (-r)-S 2j2 (-r)]. 
From Eqs. flAll) it is clear that In 5 cancels in the sum Y2j=o^j- 



(Al) 



(A2) 



(A3) 
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APPENDIX B: EXPANSIONS OF T AND V FUNCTIONS ABOUT HALF- 
INTEGER ARGUMENTS 

In this appendix, we present some useful relations between derivatives of the ip func- 
tion with half-integer arguments and the ip and f3 functions with integer arguments, and 
consider the expansion of the V function in the vicinity of half-integer arguments. 

Starting from the well-known relations between the ip and (3 functions, 




2- 



(Bl) 



and differentiating them m (m > 0) times, we have 




(B2) 



where ip^ m '(z) denotes the m-th derivative of il)(z) etc. We can combine Eqs. fIBlj) and 
(jB2p as 



^ m \z) = 2 m [^ m \2z) - ^ m \2z)] - 5 0m l 2} 
^ ( Z + l) = ^ ^ im) ^ + ^ (m) (2^)] - <Wa, (B3) 



where 5 mn is the Kronecker symbol. 



Using the series representations of the tp and (3 functions 45j, 



oo 



1 



V>(1) + (z - 1) £ 



oo 



E 



k + z' 



(B4) 



fc=0 



we obtain the following relations: 



V>(n + l)=V(l) + Si(n), 
^){n + 1) = (-l) m m\[S m+1 (n) - Cm+il 

/3(n + l) = (-l) n [Z 2 +S-i(n)], 
pW(n + 1) = (-l) m+n m![S- (m+ i)H - S- (m+ i)(oo)] 



(B5) 
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where S m (n) is defined in Eq. (1361) . 

Thus, Eqs. ( 1B3j) and (1B5j) lead to the following results for the "sums" S m with half- 
integer arguments 511 ]: 



& (|) = Si(n) + (-l) n S-i(n) - [1 - 



On 



2 m - i [^ m (n) + (-l) n S_ m (n)] + [1 - (-l) n ](l 



im-l 



(m > 2). (B6) 



These equations are useful for expansions of the F function in the vicinities of half-integer 
arguments. Indeed, using a well-known formula for the expansions of the T function about 
integer values, which was used, e.g., in Ref. [461 ] . 



T(n + 1 + 5) 
n\T(l + S) 



exp 



1 1 

oo 



X^*(n)(-*)' 



fc=i 



(B7) 



where 7^ is Euler's constant, we find the corresponding expansions about half-integer 
values to be 

00 

exp 



r(n/2 + 1 + 5) 



00 1 
fc=i 



(B8) 



r(n/2 + 1)T(1 + S) 

where S m (n/2) is given by Eq. flB6j) . Such expansions are useful in many applications, 
including those in Ref. 47| and references cited therein. 



APPENDIX C: TRANSFORMATIONS OF POLYLOGARITHMS OF 
WEIGHT FOUR 

In this appendix, we present relations between the generalized polylogarithms S a ,b of 
weight four (a + b = 4) with different arguments. Transformations at lower weights can be 



found in the literature 32|. Although the derivation of these formulas is straightforward, 
we present them here for the convenience of interested readers. At weight four, there are 
three independent Nielsen polylogarithms, which we choose to be Li 4 , Si,3, and S2,2- 
1. Relations for the functions with arguments 1 — y and y: 

Li 4 (l -y)=U- Si,s(y) + ln(l - y)[Ca - Si, 2 (y)] + \ hi 2 (l - y) [C2 - U 2 {y)\ 

-\\v?(l-y) hxy, 
o 

S 2 , 2 (l ~y) = \U - SzM + InySiM + M± ~ 2/)[Cs - Lig(j/) + InyLi 2 (y)] 
+ ^ln 2 (l-y) ln 2 y, 

Si, 3 (l - y) = C4 - LU(y) + ^ y Li 3 (y) - ^ In 2 y Li 2 (y) - J In 3 y ln(l - y). (CI) 

2 o 
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2. Relations for the functions with arguments — 1/y and — y: 
Li 4 (-^) = -L14H/) - ^4 - ^hi 2 ?/- ^ln 4 y, 

5 2 . 2 (-^) = S 2 , 2 (-y) - 2Li 4 (-y) - ^4 - lny[C 3 - Li 3 (-y)] + ^ ln 4 y, 

51.3 (-^) = ~ Si,s(-y) + S 2)2 (-y) - LU(-y) - C4 - lny[Si, 2 (-y) - Li 3 (-y)] 

-Iln 2 yLi 2 (-y)-lln 4 y. (C2) 

3. Relations for the functions with arguments (y — l)/y and y: 

Li 4 (^j^) = U M + Si, 3 (y) - S 2)2 (y) - \<U + ln(l - y)[Si, 2 (y) - Li 3 (y)] 
+ i ln 2 (l - y) Li 2 (y) - ^2 In 2 ^ + ^ ln 4 (l - y) 

-lln 4 ^, 

24 y ' 

S 2 ,2 y-^-) = 2 U M - S 2)2 (y) - ^4 + my Si, 2 (y) + In l—^G - ln[(l - y)y] Li 3 (y) 

V y j y 

+ In y ln(l - y) Li 2 (y) + J In 2 y ln 2 (l - y) - i In 3 y ln(l - y) 

2/-^ _t; a ,„..t: /.a , 1 ,^2,. T: /„ ,\ , 1,3, 



Si,3 y-^T) = Uii3j) ~ C4 ~ ln y Lis(y) + 2 y Li2(2/) + 6 y ln(1 ~ y) 

-^In 4 !/- (C3) 
4. Relations for the functions with arguments y/(y — 1) and y: 

Li 4 (^-y) = S 2 , 2 (y) - Li 4 (y) - Si, 3 (y) + ln(l - y)[Li 3 (y) - Si,2 
-iln 2 (l-y)Li 2 (y)-lln 4 (l-y), 



S2 < 2 = S2 ' 2 (^) " 2 Sl -3(f) - ln ( X " + 2i ln4(1 " f )' 

Si,3 (^l)=-Si,3(y)-^ln 4 (l-y). 



(C4) 
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5. Relations for the functions with arguments 1/(1 + y) and — y: 
1 \ . , s 1 



'1,3 



Li4 ( T+^J = Sl - 3(-2/) + C4 + ln(1 + ^)[ Sl - 2 (-^) " Ca] + g ln ( X + ?/)[ Li 2(-?/) + Ca] 
+ ^ln 3 (l + y) lny-^ln 4 (l + y), 

S2 - 2 ( i^) = 2 M"?/) " S 2 , 2 (-y) + \u + + ln[y(l + j/)] Si, 2 (-3/) 

- ln(l + y) [Li 3 (-y) + Cs] + hi(l + y) lny Li 2 (-y) + ~ ln 2 (l + y) In 2 y 

- iln 3 (l + y) lnj/ + ^ln 4 (l + ?/), 



1 + 3/ 



Si, 3 (-y) - S 2|2 (-y) + Li 4 (-y) + C 4 + lny[Si, 2 (-y) - Li 3 (-y)] 
+ \ In 2 y Li 2 (-j/) + 1 In 4 y - 1 In 4 1±^. (C5) 



Equations (ICip and (IC2p were directly obtained from Ref . [32j , where they are presented 
for the generalized polylogarithms S a ,6 with arbitrary values of a and b, but in some 
complicated form less convenient for applications. Equations (1C3j) - (1C5j) were found by 
iterated application of Eqs. (IClj) and (1C2I) and equations from Ref. [32|. They are sim- 
ple and useful for applications together with equations for S a ,b from Ref. [32J, with the 
constraints a + b = 2 or a + b = 3. 
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